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We numerically investigate the influence of interactions on the generalized Hong-Ou-Mandel 
(HOM) effect for bosonic particles and show results for the cases of = 2, = 3 and A'' = 4 

bosons interacting with a beam splitter, whose role is played by a 5-barrier. In particular, we focus 
on the effect of attractive interactions and compare the results with the repulsive case, as well as 
with the analytically available results for the non-interacting case (that we use as a benchmark). 

We observe a fermionization effect both for growing repulsive and attractive interactions, i.e., the 
dip in the HOM coincidence count is progressively smeared out, for increasing interaction strengths. 

The role of input asymmetries is also explored. 
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I. INTRODUCTION 

When two indistinguishable photons simultaneously 
hit a 50 % - 50 % beam splitter, one in each input port, 
they interfere destructively and no coincidence counts of 
particles detected in both output ports can be observed. 
This inherently quantum mechanical effect is a striking 
manifestation of the bosonic properties of the photons. 
Since its first observation in 1987, the celebrated Hong- 
Ou-Mandel (HOM) [I| effect has triggered a multitude of 
investigations and has been the focus of ongoing research: 
it has found application in the creation of post-selected 
entanglement between photon pairs 0, and can be ex¬ 
ploited for logic gates in linear optical quantum compu¬ 
tation Q, or to demonstrate the purity of a solid-state 
single-photon source Q. The HOM effect has also been 
experimentally demonstrated with single photons emit¬ 
ted by two independently trapped single atoms The 
influence of varying properties of the beam-splitter on 
the joint probability distribution (of finding both photons 
on the same side) has been investigated in Q. Experi¬ 
mental realizations have been extended to larger particle 
numbers, namely three photons impinging on a multi- 
port mixer Q, and to one and two-photon pairs Q. A 
review on multi-photon experiments and the generalized 
Hong-Ou-Mandel effect is given in Q. 

More recently, the HOM effect has been generalized to 
massive particles, i. e. to N bosons or fermions passing 
simultaneously through a symmetric Bell multiport beam 
splitter [T^ , and to a large number of particles impinging 
on a single beam splitter 0[i3- In the latter case it has 
been shown for the balanced beam-splitter that if an even 
number of particles impinges on the beam-splitter from 
either side, an even number must also emerge from each 
side. A recent proposal suggests to observe the HOM 
effect with colliding Bose-Einstein condensates [l^ for 
a set-up that has already been used to demonstrate the 
violation of the Cauchy-Schwarz inequality [l^ . 


Effects of the interparticle interaction have been in¬ 
vestigated in [I^ within a Bose-Hubbard set-up where 
a transition from bunching to antibunching can be ob¬ 
served for growing repulsive interparticle interactions, 
corresponding to a fermionization of the particles [il[i3- 
Bosonic atoms have also been proposed as a suitable can¬ 
didate for Knill-Laflamme-Milburn quantum computa¬ 
tion, with the advantage of very controllable input state 
preparation [I^. It is relevant to highlight here that 
two-particle quantum interference akin to the Hong-Ou- 
Mandel effect has recently been observed in the context of 
independently prepared bosonic atoms in tunnel-coupled 
optical tweezers [l9| . The latter experiment follows up on 
another experimental realization of such quan tum inter¬ 
ference effects in the context of electrons [^. It should 
thus be clear that the remarkable advances in the set ting 
of ultracold gases enable the very accurate counting 
and controllability [12, [12| of bosonic atoms towards such 
experiments. Motivated by the fact that the parity of the 
number of atoms in a potential well is experimentally 
accessible more straightforwardly than the exact parti¬ 
cle number it has been pointed out that parity measure¬ 
ments can yield useful signatures of the generalized HOM 
effect [l2| ■ Recent experimental results even demonstrate 
the simultaneous determination of the number of atoms 
in each well of a double-well trap with single-atom reso¬ 
lution for up to AI = 500 atoms per well [^ . 

Here, we will chiefly focus on the effect of attractive 
interparticle interactions on the HOM dip and demon¬ 
strate that also in this case a transition from bunching 
to antibunching with growing interaction strength can 
be observed. This can again be interpreted as a fermion¬ 
ization effect, in agreement with [25l| . where the Bose- 
Fermi mapping has been demonstrated for attractive ID 
bosons. The latter generalization enables the formation 
of gas-like states that fermionize as the attraction be¬ 
comes stronger. It is that, somewhat counter-intuitive 
(on the basis of the nature of the interaction) feature that 
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we also find in the present setting. Our investigation fo¬ 
cuses on the (generalized) HOM effect for a ID Bose gas 
on the A^-particle level, enabling true quantum behavior. 
We focus on the cases N = 2, N = 3 and N = A. In ad¬ 
dition to examining the somewhat less studied attractive 
case, we compare the results with those of the repulsive 
case and importantly with the non-interacting case that 
can be addressed, in principle, in its full generality for 
arbitrary iV, as will be discussed below. 

In the original HOM experiment the photons emerge 
as a result of quantum interference in a superposition of 
states |2,0) and |0, 2) and thus in a measurement would 
always be found on the same side. As an aside, we note in 
passing, effectively regarding the large N limit of mean- 
field matter waves that a recent proposal suggests an ana¬ 
logue of the HOM effect with bright solitons [2^ . In this 
classical case it is found that the indistinguishability of 
the particles yields a 0.5 split mass on either side for soli¬ 
tary waves (by parity symmetry at this mean-field level). 
But for very slight deviations it can be observed that 
all the particles always end up on the same side of the 
barrier potential. Recently, also the collisional dynamics 
of matter-wave solitons (in the absence of a barrier po¬ 
tential) have been investigated experimentally [l^. Ad¬ 
ditionally, the setting of interactions of individual sin¬ 
gle [ 2 ^ and multi-component [ 2 ^ solitary waves with 
barriers (that play the role of the beam splitter here) is 
certainly within experimental reach. Collisions have also 
been suggested as a way to create entanglement between 
indistinguishable solitons and initially independent 
and indistinguishable boson pairs [3l|. Collisions of dis¬ 
tinguishable solitons have been proposed as a possibility 
to create mesoscopic Bell states HI. 

Our presentation is structured as follows: Sec. El gives 
an introduction to the model system, the Lieb-Liniger(- 
McGuire) model, and a quantum beam-splitter. In 
Sec.EIl we present numerical results for the collision of 
two monomers, two dimers and the asymmetric case of a 
monomer and a dimer on an additional barrier potential. 
Sec.HV] summarizes our findings and presents a number 
of future challenges. 


II. MODEL SYSTEM 

The iV-particle dynamics of interacting bosons in 
(quasi-)one-dimensional geometries (corresponding to a 
3D geometry with tightly confined radial degrees of free¬ 
dom) can be described within the exactly solvable Lieb- 
Liniger (-McGuire) model [33,[33 

N ,2 ^-1 ^ 

Hll = + 9S{xj-Xn). 

i=l i=l n=j+l 

( 1 ) 

Here, contact interaction between the N bosons is as¬ 
sumed and quantihed with the interaction parameter g. 


In the following we will investigate the scattering dynam¬ 
ics within an additional harmonic confinement (emulat¬ 
ing the typical parabolic trap relevant to experimental 
settings; cf. [^[2^) at a repulsive delta-like barrier po¬ 
tential in the middle of the harmonic confinement, yield¬ 
ing the full Hamiltonian 


N 

i? = i?LL+5Z^-t(x,) (2) 

i=i 

with the external potential Vext{xj) = +voS{xj). 

For the numerical implementation the system can be 
discretized via the Bose-Hubbard Hamiltonian 

.^discretized = ~ J -(- — ^ ^ hj {flj — 1 ) 
3 j 

+A njj^ + voSjfi, (3) 

j 

where denotes the annihilation (creation) operator 

for lattice site j, fij = cjcj is the particle number oper¬ 
ator, and U the on-site interaction strength. The tun¬ 
neling strength is given by J ~ h^/2mb^ with lattice 
spacing 5 for 5 —>■ 0, A = defines the strength 

of the harmonic confinement and vq the strength of the 
delta-like barrier potential. For sufficiently small lattice 
spacing 5 —>■ 0 the Lieb-Liniger model (HI) with additional 
harmonic confinement is recovered. 


A. Quantum beam-splitter 


The repulsive delta-like barrier potential acts as a 
beam-splitter. In this section we present a theoretical 
description for the noninteracting case, following 
(see also [HI, [l3)- Assume a generalized beam-splitter 
with two incoming modes (oi) and ( 02 ) being coupled at 
the beam-splitting device and possible additional phase 
shifts, as depicted in Fig. [TJ This system is described by 
the time-evolution operator 


Uc(0,^) 


/ cos(0/2) ie®‘^ sin(0/2)\ 

Vfe-*'^sin(6l/2) cos(6»/2) J ’ 


( 4 ) 


where 0 is related to the complex transmission and reflec¬ 
tion coefficients t = cos(0/2) and r = jsin(0/2), obeying 
-I- |rp = 1 and tr* + t*r = 0. It will be seen that the 
occupation probabilities of the possible output states are 
independent of the phase cj). 

We first review the well-known case of a balanced 
beam-splitter with one particle in each input mode. The 
output state then reads 


C/c(7r/2,0)|l,l) (5) 

(aj) -k a\al - 44 + * (4) ^ |0) 0) (6) 
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For the bosonic case the application of the corresponding 
commutation relations, [ 0 ^, 0 ]] = = Sij and 

[a|^,aj] = [a^,a^] = 0, yields the charateristic bunching: 

t/e(V2,0)|l,l) = ^(|2,0) + |0,2)). (7) 

Both particles would always be detected in the same 
output mode. For fermions obeying the commutation 
relations { 0 ^, 0 ]} = a^a] + a]aj = % and {a-,a]} = 
{ 0 ^, 0 ^} = 0, the output state exhibits anti-bunching: 

C/c(V2,0)|l,l) = -|l,l). (8) 

For bosons this calculation can be straightforwardly ex¬ 
tended to larger particle numbers: consider a Fock input 
state with n particles in mode (oi) and m particles in 
mode ( 02 ). The output state then is given by the super¬ 
position state 

Uc{0,(())\n,m) 

/ n\ fkfra — k + n — Vs 

I ; V rn-k ) 

X (z sin(6l/2))™-'=+' cos(0/2)'^+"-' exp(t^(m - k - 1)) 
x\m — k + n — Vk + 1). (9) 

This corresponds to the occupation probability 


P\nf,mf) 

= \{nf,rnf\Uc{9,(j))\n,m)\'^ 

= n\m\nf\mf\ cos(6»/2)2(™/+") sin(6»/2)2(™-'"/) 

(—1)* sin(6l/2)^^ cos(0/2)“^^ 


E' 

/=0 


l\(n — l)\{mf — l)\(m — nif -|- /)! 


,( 10 ) 


of a final state in agreement with Eq. (14) 

from m for the balanced beam-splitter. A derivation 
of the occupation probabilities for arbitrary reflectiv¬ 
ity/transmittivity has also been given in [T^ . The prime 
indicates that summands resulting in undefined negative 
factorials do not contribute. 

For indistinguishable fermions on the other hand an 
extension of this description to iV > 2 is not directly 
possible for this two-mode set-up due to the Pauli prin¬ 
ciple: an occupation of any mode with more than one 
particle is not allowed. Note that also for the description 
of N fermions in a multi-mode setting as m the number 
of input / output modes has to be equal to or larger than 
the number of particles. 


III. NUMERICAL RESULTS 

Our numerical results are obtained by implement¬ 
ing ([ 3 ]) for sufficiently small lattice spacing b, to emulate 
the Lieb-Liniger model. First, we investigate the direct 
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FIG. 1: Schematic depiction of a quantum beam-splitter (BS) 
with input modes (ai) and ( 02 ), detectors (D) at both outputs 
and optional additional phase shifts. 





vomXiJh^ \U/J\ 


FIG. 2: (Golor online) Analogue of Hong-Ou-Mandel effect for 
two monomers colliding on a narrow delta-like potential bar¬ 
rier situated in the middle of an additional harmonic conhne- 
ment with oj = 1.0. The monomers are initially displaced by 
equal amounts |2;o| to the left and right of the barrier poten¬ 
tial. Shown is the occupation probability P|i,i) of state |1,1) 
at t/T = 0.5 vs. the scaled height vomX^^^^fh^ of the bar¬ 
rier potential for interaction strengths U/ J = 0.0 (solid black 
line), Uj J = —0.5 (solid blue line) and U/ J = —1.0 (solid red 
line) and different initial displacements (a) |a:o|/Aoac.i = 5, 
(b) |a:o|/Aosc,i = 10 and (c) |xo|/Aosc,i = 15 of the monomers. 
Results for repulsive interactions U j J = 0.5 (dashed blue 
line) and U/J = 1.0 (dashed red line) are displayed as well. 
Dashed green line: analytical prediction given by Eq. (nsj. 
(d) Occupation probability of state |1,1) vs. ratio |U/J| of 
interaction strength to tunneling strength for attractive (black 
line) and repulsive interactions (red line) for initial displace¬ 
ments Isol/Aoac.i = 10. 
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FIG. 3: (Color online) Influence of shifts a;o,sh/Aosci around 
the initial position lo/Aoaci = 10 of the right particle for 
monomers on the occupation probability of state |1,1). The 
height = 200 is chosen to ensure 50% — 50%- 

splitting in the case of no shift. Black: U/J = 0.0, blue: 
U/J = —1.0 (solid line), U/J — 1.0 (dashed line), red: U/J = 
— 2.0 (solid line), U/J = 2.0 (dashed line). Same parameters 
as in Fig. [21(b). 

analogue to the original HOM experiment [if. We simu¬ 
late the collision of two single particles on a barrier po¬ 
tential situated in the middle a; = 0 of an additional har¬ 
monic confinement. The initial state is given as the (sym¬ 
metrized) product state of two single particles, each of 
them initially prepared in well-separated harmonic con¬ 
finements with same trapping frequency and potential 
minima located at a; = Txq. Subsequently, the harmonic 
confinements are replaced by a single harmonic conhne- 
ment with potential minimum at a; = 0 implying that 
both particles then are located at the classical turning 
points of the harmonic trap and start to accelerate with 
oscillation period T = ^ towards the trap center where 
they simultaneously hit the narrow, delta-like potential 
barrier. In Fig. [2] we display the occupation probability 
of state 11,1) at t = T/2 for different interaction and 
potential strengths and different initial displacements re¬ 
sulting in different kinetic energies when the particles 
hit the barrier at t = T/A. Without interparticle inter¬ 
actions we observe a behavior analogous to the original 
Hong-Ou-Mandel experiment: if the potential height is 
chosen to ensure 50%-50%-splitting for each single par¬ 
ticle, we clearly observe the HOM dip with a negligi¬ 
ble occupation probability of state |1,1); the system is 
found in the superposition state (I^), nicely illustrating 
the bosonic symmetry and interference properties of the 
particles, in the absence of interactions. Again analogous 
to the original HOM effect changes in potential height, 
resulting in deviations from 50%-50%-splitting, lead to a 
growing contribution of |1,1), in agreement with the an- 


FIG. 4: (Color online) Same as Fig. [2] for two dimers and 
uj = 4.0. Initial displacement by 2.2Aosci • Displayed is the 
occupation probability of the states |2,2), |3,1) and |4, 0) at 
t/T = 0.5 vs. the scaled height uomAosc^/fi^ of the bar¬ 
rier potential for interactions strengths U/J = 0.0 (black), 
U/J = —0.5 (blue), U/J — —1.0 (red), U/J — 0.5 (dashed 
blue) and U/J = 1.0 (dashed red). Vertical dotted black line: 
Value of potential height for 50%-50%-splitting. Horizontal 
dotted black line: Analytical expectation for the linear case 
and a balanced beam splitter. Note that P| 3 ,i) = P| 3 ,i) = 0 
is the analytical prediction for this case. Dashed green line: 
analytical prediction given by Eq. m- The slight deviation 
from the analytical prediction is a remaining lattice effect due 
to the restriction of the numerics for A = 4 to a smaller num¬ 
ber of lattice sites, here 75. 


alytical prediction (uni). Notice that both for this and for 
all the cases that will follow, the difference between the 
numerically obtained solid black line and the analytically 
derived green dash-dotted line on the basis of Eq. (nni) 
will be a measure of the quality of our approximation 
(via the Bose-Hubbard underlying lattice model) of the 
continuum setting. In all the case examples displayed 
below, these deviations are very small. 

For repulsive interactions a transition from bunching 
to antibunching has been predicted with growing inter¬ 
actions [l^ . A particularly intriguing feature of our find¬ 
ings is that also in the presence of attractive interactions 
the Hong-Ou-Mandel dip is suppressed and a transition 
from bunching to antibunching can be observed. Such a 
somewhat counter-intuitive (on the basis of the attrac¬ 
tive interactions) phenomenon is reminiscent of similar 
fermionization features as those reported in [2^; cf. Figs, 
like Fig. Ib therein. For comparison, Fig.|2]shows results 
both for attractive and for repulsive interactions. It can 
be seen that for increasing kinetic energy of the particles 
the suppression of the HOM dip is reduced, cf. Fig.|2](b) 
and Fig. m{c). This is natural to expect since in the latter 
setting, the kinetic effects dominate the interaction ones 











5 


and progressively (as the kinetic energy increases) restore 
the non-interacting picture. While differences between 
repulsive and attractive interactions are found for lower 
kinetic energies attractive and repulsive interactions yield 
quantitatively similar results for larger kinetic energies. 
This can be explained by the shorter interaction time in 
the middle of the harmonic trap for larger kinetic ener¬ 
gies, rendering interactions less influential. Summariz¬ 
ing these results, the HOM dip is suppressed with grow¬ 
ing attractive interaction strengths for all cases, but this 
suppression is most effective for low center-of-mass ki¬ 
netic energies. Fig. [5] (d) compares the suppression of 
the HOM dip for continously growing attractive and re¬ 
pulsive interaction strengths. It can be seen that up to 
|17/J| ~ 1 the occupation probability of state |1,1) is the 
same but for larger interaction strengths the HOM dip is 
suppressed equally or in some cases slightly more strongly 
in the presence of attractive interactions (in comparison 
to repulsive ones). For repulsive interactions the |1,1) 
configuration should be energetically favored, while for 
attractive interactions the superposition state ([7]) should 
intuitively be favored. We attribute our observation of 
the antibunching to the fact that, in analo gy t o the Bose- 
Fermi mapping in the repulsive case [l^ . [T^ (for large 
interactions U), also attractive bosons in ID undergo a 
fermionization for growing interaction strengths [^ . 

At the mean-field level, applicable for large N, a 0.5 
split mass on either side is the result for a balanced beam¬ 
splitter (due to parity symmetry). Yet, it has recently 
been shown that for very slight deviations from indis- 
tinguishability, it can be observed that almost all the 
particles may end up on the same side of the barrier po- 
tential, constituting a mean-field analogue of the HOM- 
effect [26|. This observation prompted us to investigate 
the role of asymmetries in the set-up by introducing 
slightly different initial displacements of the monomers, 
cf. Fig. [3] due to the harmonic confinement the parti¬ 
cles then still hit the barrier at the same point of time 
but with different kinetic energies. While on the mean- 
field level such asymmetries have caused the occurrence 
of HOM-like behavior, here we observe a complementary 
behavior: at the few-particle quantum-mechanical level 
such asymmetries result in a suppression of the HOM 
dip. For very large interaction strengths \U/J\ « 2 the 
occupation ppp) remains close to its already large value 
of 0.51 (repulsive) and 0.55 (attractive) without shifts of 
the initial position. This suggests that as N is increased 
a transition should arise between the small-behavior 
presenting a dip at vanishing as ym metry for N = 2 and 
odd-even oscillations for > 2 and the large N be¬ 
havior bearing a 0.5 splitting fraction at vanishing asym¬ 
metry and a potentially close to unity fraction for slight 
asymmetries. Unfortunately, this regime is not currently 
accessible to our techniques, yet would constitute an in¬ 
triguing theoretical (and experimental) problem to con¬ 
sider in future work. 

We now turn to the investigation of the situation of two 
dimers, cf. Fig.Ul In this case, the ground state of each of 



FIG. 5: (Color online) Same as Figs. [2] and [4] and for a dimer 
(initially left) and a monomer (initially right) and ta = 1.5. 
Initial displacement by 3.3Aosci. Displayed is the occupa¬ 
tion probability of the states |3,0), |2,1), |1,2) and |0,3) 
at t/T = 0.5 vs. the scaled height of the bar¬ 

rier potential for interactions strengths U/J = 0.0 (black), 
U/J = —0.5 (blue), UjJ — —1.0 (red), UjJ — 0.5 (dashed 
blue) and UjJ = 1.0 (dashed red). Calculated on 121 lattice 
points. Vertical dotted black line: Value of potential height 
for 50%-50%-splitting. Horizontal dotted black line: Ana¬ 
lytical expectation for the linear case and a balanced beam 
splitter. Dashed green line: analytical prediction given by 
Eq. (IIOII . The slight deviation from the analytical prediction 
again is a residual lattice effect. 


the dimers is obtained by imaginary time evolution: the 
many-particle initial wave function is taken as a prod¬ 
uct wave function of both dimer ground states. While 
for the monomers a large number of lattice sites was ac¬ 
cessible, for the dimers the numerically accessible size of 
the Hilbert space is considerably reduced. We choose 75 
lattice sites. Once again, we utilize the non-interacting 
case as a benchmark for our results, finding good agree¬ 
ment with the analytical prediction of Eq. (fTOll . The 
slight deviation is caused by residual lattice effects. Due 
to symmetry we only display occupation probabilities of 
the states |2,2), |3,1) and |4,0). It can be seen that for 
the balanced beam-splitter only states with even particle 
numbers on either side of the beam-splitter contribute to 
the output state. This corresponds to the predicted odd- 
even oscillations for the generalized HOM effect in [HI- 
Fig. |4] also displays the effect of interactions. It can be 
observed that the Fock-states are affected in a different 
manner: while the occupation probability of state |4, 0) 
(and |0,4)) is affected in a similar way, regardless of the 
sign of the interaction, considerable differences can be 
observed for the other states. Importantly, we note that 
in this case no direct analog of the Hong-Ou-Mandel ef¬ 
fect exists, as even in the absence of interactions and for a 
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50 % - 50 % beam-splitter, the |2,2) state is not fully sup¬ 
pressed (although the asymmetric |3,1) and |1,3) ones 
are, as indicated above). What our results suggest is 
that as we depart from this limit, the preferred non¬ 
interacting state of |4,0) gets suppressed in favor of in¬ 
creased probabilities of both the |2, 2) and |3,1) states. 
The attractive case progressively leads to fermionization 
in the form of the |2, 2) state, while repulsive interactions 
progressively lead to an increased probability of the |3,1) 
state. 

As discussed in Sec. Ill Al the analytical description 
in terms of the two-mode quantum beam-splitter for 
fermions cannot be extended to cases with N > 2. While 
this prevents a direct comparison of our results with the 
fermionic case, we nonetheless presume that the effects 
of interaction we observe can be attributed to an effec¬ 
tive fermionization: Ref. [2^ demonstrates the occur¬ 
rence of fermionization for an attractive one-dimensional 
Bose gases by the formulation of a Bose-Fermi mapping 
and discusses the stationary states of a few bosons in a 
harmonic confinement. 

Finally, we investigate the asymmetric case of = 3, 
the collision of a dimer and a monomer, cf. Fig. where 
again we interpret the result of interactions as suggest¬ 
ing an effective fermionization, following a similar rea¬ 
soning as above. As the results for the interacting case 
are governed by a directional effect, we display the oc¬ 
cupation probabilities of all possible final states. The 
non-interacting case again is in accordance with the ana¬ 
lytical prediction (fTOl) . In the absence of a barrier, here, 
we get solely the |1, 2) state, while under a 50 % - 50 % 
splitting (in the non-interacting case) |3, 0) and |0, 3) are 
equally populated and favored in comparison to |1, 2) and 
|2,1) (who are also equally probably between them). The 
interactions generally favor less the bunched states |3, 0) 
and |0, 3), although the different interaction signs lead to 
a more pronounced effect in the former in comparison to 
the latter. Similarly, the interactions lead to an increased 
probability of |1, 2) and |2,1), although again the differ¬ 
ence in signs affects more the former in comparison to 
the latter. 


IV. CONCLUSION 

We have numerically investigated the generalized 
HOM-effect for bosonic atoms in a (quasi-)one- 
dimensional geometry with a special emphasis on the 
influence of interactions. In accordance to existing litera¬ 
ture for a repulsive Bose-Hubbard set-up [T^ , we observe 
a transition from bunching to antibunching for repulsive 
interactions. In addition, we have presented numerical 
results that show such a fermionization and generally 


similar anti-bunching trends also for attractive interac¬ 
tions. Our results have been illustrated beyond the well- 
known case of V = 2 atoms, for those with N = 3 and 
N = A. The role of the strength of the interactions (and 
its relative influence in comparison to the kinetic energy) 
was illustrated and the influence of a potential asymme¬ 
try in the initial configuration was also discussed. We be¬ 
lieve that this presentation, thus, yields a fairly complete 
picture of the generalized form of HOM-type (numerical) 
experiments for small atom numbers N. 

The high controllability available in current experi¬ 
ments with atoms, as discussed in the Introduction, ap¬ 
pears to render very plausible and accessible an experi¬ 
mental observation of the generalized Hong-Ou-Mandel 
effect, exploring, e.g., odd-even oscillations in the occu¬ 
pation of the final state in the non-interacting limit, as 
well as suppression of HOM-like dips for increasing inter¬ 
action strengths, as discussed herein. The control of in¬ 
teraction strengths via tools such as magnetically or 
optically [37| induced Feshbach resonances should enable 
the systematic observation of the numerically obtained 
features. We argue that for larger displacements out of 
the trap center the interactions play a negligible role and 
results close to the linear case can be obtained. Due to 
their stability, bright solitons are a particularly promis¬ 
ing candidate to observe the odd-even oscillations in the 
regime of low particle numbers where quantum effects 
are important. Additionally, a systematic exploration of 
regimes of progressively larger V, to explore the system- 
atics of the transition from the small N to the mean-field 
behavior would be of particular interest for future stud¬ 
ies. 
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